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Abstract The properties of hairy black holes in Einstein- Yang-Mills (EYM) theory 
are reviewed, focusing on spherically symmetric solutions. In particular, in asymp- 
totically anti-de Sitter space (adS) stable black hole hair is known to exist for su(2) 
EYM. We review recent work in which it is shown that stable hair also exists in 
su(A^) EYM for arbitrary A^, so that there is no upper limit on how much stable hair 
a black hole in adS can possess. 



1 Introduction 

We begin by very briefly reviewing the "no-hair" conjecture and motivating the 
study of hairy black holes. 



1.1 The "no-hair" conjecture 

The black hole "no-hair" conjecture II141I states that (see, for example, 1151115211771 
FTSI [791 ITTtI for detailed reviews and comprehensive lists of references): 



All stationary, asymptotically flat, four-dimensional black hole equilibrium 
solutions of the Einstein equations in vacuum or with an electromagnetic field 
are characterized by their mass, angular momentum, and (electric or mag- 
netic) charge. 
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According to the no-hair conjecture, black holes are therefore extraordinarily 
simple objects, whose geometry (exterior to the event horizon) is a member of the 
Kerr-Newman family and completely determined by just three quantities (mass, 
angular momentum and charge). Furthermore, these quantities are global charges 
which can (at least in principle) be measured at infinity, far from the black hole 
event horizon. If a black hole is formed by the gravitational collapse of a dying star, 
the initial star will be a highly complex object described by many different param- 
eters. The final, equilibrium, black hole is, by contrast, rather simple and described 
by a very small number of quantities. During the process of the formation of a black 
hole, an enormous amount of (classical) information about the star which collapsed 
has therefore been lost. Similarly, if a complicated object is thrown down a black 
hole event horizon, once the system settles down, the only changes in the final state 
will be changes in the total mass, total angular momentum and total charge. Ad- 
vances in astrometry 01731 and future gravitational wave detectors Q may even be 
able to probe the validity of the "no-hair" conjecture for astrophysical black holes 
by verifying that the mass, angular momentum and quadrupole moment Q2 of the 
black hole satisfy the relation Q2 ~ J^/M which holds for Kerr black holes. 

The "no-hair" conjecture, stated above, has been proved by means of much 
complicated and beautiful mathematics (as reviewed in, for example, lISTl l52l ITTl 
l78l |79l II 171 ). subject to the assumptions of stationarity, asymptotic flatness, four- 
dimensional space-time and the electrovac Einstein equations. It is perhaps unsur- 
prising that if one or more of these assumptions is relaxed, then the conjecture does 
not necessarily hold. For example, if a negative cosmological constant is included, 
so that the space-time is no longer asymptotically flat but instead approaches anti- 
de Sitter (adS) space at infinity, then the event horizon of the black hole is not 
necessarily spherical, giving rise to "topological" black holes (see, for example, 
lfT8ll64ll971II02l II III 11641 ). More recently, the discovery of "black ring" solutions 
in five space-time dimensions ( II6OI . see IMI for a recent review) and the even more 
complicated "black Saturn" 1591 solutions indicates that Einstein-Maxwell theory 
has a rich space of black solutions in higher dimensions, which are not given in 
terms of the Myers-Perry 11201 metric (which is the generalization of the Kerr- 
Newman geometry to higher dimensions). 



1.2 Hairy black holes 

In this article we consider what happens when the other condition in the "no-hair" 
conjecture, namely that the Einstein equations involve electrovac matter only, is re- 
laxed. The "generalized" version of the no-hair conjecture ||79l states that all station- 
ary black hole solutions of the Einstein equations with any type of self-gravitating 
matter field are determined uniquely by their mass, angular momentum and a set of 
global charges. Even in asymptotically flat space, this conjecture does not hold, even 
for the simplest type of self-gravitating matter, a scalar field. The first such counter- 
example is the famous BBMB black hole IIT2I [131 l27l which has the same metric 
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as the extremal Reissner-Nordstrom black hole but possesses a conformally cou- 
pled scalar field. However, this solution is controversial due to the divergence of the 
scalar field on the event horizon ||157| and is also highly unstable ||48|. Therefore, 
in some ways the first "hairy" black hole is considered to be the Gibbons solution 
ITTTI . which describes a Reissner-Nordstrom black hole with a non-trivial dilaton 
field. While there are many results which rule out scalar field hair in quite general 
models, particularly in asymptotically flat space-times (see, for example, lfT4l for a 
review), in recent years many other examples of black holes with non-trivial scalar 
field hair have been found. For example, minimally coupled scalar field hair has 
been found when the cosmological constant is positive 11601 or negative II161I and 
non-minimally coupled scalar field hair has also been considered (see, for example, 
II175II1761 and references therein). 

In this short review, we will focus on another particular matter model, Einstein- 
Yang-Mills theory (EYM), where the matter is described by a non-Abelian (Yang- 
Mills) gauge field. It is now well-known that this theory possesses "hairy" black 
hole solutions, whose metric is not a member of the Kerr-Newman family (see 
II 1701 for a detailed review). Furthermore, unlike the Kerr-Newman black holes, 
the geometry exterior to the event horizon is not determined uniquely by global 
charges measureable at infinity, although only a small number of parameters are 
required in order to describe the metric and matter field (see section [3] for further 
details). All of the asymptotically flat black hole solutions of pure EYM theory 
discovered to date are unstable fiTl (however, there are examples of asymptotically 
flat, stable hairy black holes in variants of the EYM action, such as Einstein-Skyrme 
ElllsSlllOllSlI, Einstein-non-AbeHan-Proca I73l[l09l[l58l[l59l[l62l and Einstein- 
Yang-Mills-Higgs |[T] theories). This means that, while the "letter" of the no-hair 
theorem is violated in this case (as there exist solutions which are not described 
by the Kerr-Newman metric), its "spirit" is intact, as stable equilibrium black holes 
remain simple objects, described by a few parameters if not exactly of the Kerr- 
Newman form (see II2TI for a related discussion along these lines). 

The situation is radically different if one considers EYM solutions in asymptoti- 
cally adS space, rather than asymptotically flat space. For su(2) EYM, at least some 
black hole solutions with hair are stable 1251 l26l [17411 . These stable black holes 
require one new parameter (see section|4]i to completely describe the geometry ex- 
terior to their event horizons. Therefore, one might still argue that the true "spirit" 
of the "no-hair" conjecture remains intact, that stable equilibrium black holes are 
comparatively simple objects, described by just a few parameters. 

One is therefore led to a natural question: are there hairy black hole solutions in 
adS which require an infinite number of parameters to fully describe the geometry 
and matter exterior to the event horizon? In other words, is there a limit to how 
much hair a black hole in adS can be given? This is the question we will be seeking 
to address in this article. 
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1.3 Scope of this article 

The subject of hairy black holes in EYM theory and its variants is very active, with 
many new solutions appearing each year. The review II170L written in 1998, is very 
detailed and thorough, and contains a comprehensive list of references to solutions 
known at that time. We have therefore not sought to be complete in our references 
prior to that date, and have, instead, chosen to highlight a few solutions (the selection 
being undoubtedly personal). Even considering just work after 1998, we have been 
unable to do justice to the huge body of work in this area (for example, the seminal 
paper Q has 172 arXiv citations between 1999 and the time of writing), and have 
instead chosen some examples of solutions. As well as II170I , reviews of various 
aspects of solitons and black holes in EYM can be found in ETl l66l l72l 11511 11521 



The outline of this article is as follows. In section|2]we will outline 5u(A^) EYM 
theory, including our ansatz for the gauge field and the form of the field equations. 
We will then, in section |3] briefly review some of the properties of the well-known 
asymptotically flat solutions of this theory. Our main focus in this article are asymp- 
totically adS black holes, and we begin our discussion of these in section |4]by re- 
viewing the key features of the su(2) EYM black holes in adS, before moving on to 
describe very recent work on su{N), asymptotically adS, EYM black holes in sec- 
tion|5] Our conclusions are presented in section|6l Throughout this article the metric 
has signature (—,+,+,+) and we use units in which 47rG = c = 1. 



2 su(A^) Einstein- Yang-Mills theory 

In this section we gather together all the formalism and field equations we shall 
require for our later study of black hole solutions. 



2.1 Ansatz, field equations and boundary conditions 

In this article we shall be interested in four-dimensional su(A^) EYM theory with a 
cosmological constant, described by the following action, given in suitable units: 



where R is the Ricci scalar of the geometry and A the cosmological constant. Here 
we have chosen the simplest type of EYM-like theory, many variants have been 
studied in the literature (see, for example, 117011 for a selection of examples). 
Varying the action ([T]i gives the field equations 



[1651 . 




(1) 
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O^D^Fy^'^V^Fyf+[A^,Fy^']■, (2) 
where the YM stress-energy tensor is 

T^v = T^rF^^Fy"- - y^vTTF;,,F^^. (3) 

In this article we consider only static, spherically symmetric black hole geometries, 
with metric given, in standard Schwarzschild-like co-ordinates, as 

ds^ = -HS^dt^ + H-^dr^ + r^de^ + r^ sin^ ed^^, (4) 

where the metric functions /i and S depend on the radial co-ordinate r only. In the 
presence of a negative cosmological constant A < 0, we write the metric function jj. 
as 

(5) 

r 3 

The most general, spherically symmetric, ansatz for the su(A^) gauge potential is 



A = £/dt + ^dr+-{C-C") de-^[(C + C")sm9+Dcos9] d(j), (6) 

where £/, SS, C and D are all (A^ x A^) matrices and is the Hermitian conjugate 
of C. The matrices ^ and S§ are purely imaginary, diagonal, traceless and depend 
only on the radial co-ordinate r. The matrix C is upper-triangular, with non-zero 
entries only immediately above the diagonal: 

C,-,,+i=a),,-(ryi'^W, (7) 

for y = 1 , . . . , — 1 . In addition, D is a constant matrix: 

£) = Diag(A'-l,A^-3,...,-A' + 3,-A'+l). (8) 

Here we are primarily interested only in purely magnetic solutions, so we set = Q. 
We may also take ^ = by a choice of gauge [98] . From now on we will assume 
that all the (Oj{r) are non-zero (see, for example, Il69ll94l l95 96] for the possibilities 
in asymptotically flat space if this assumption does not hold). In this case one of the 
Yang-Mills equations becomes 



Yj^Q V; = 1,...,A^-1. (9) 
Our ansatz for the Yang-Mills potential therefore reduces to 

A = i (C-C^) de-^- [(C + C^) sin0+Dcos0] (10) 
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where the only non-zero entries of the matrix C are 

Cj.j+i^COj{r). (11) 

The gauge field is therefore described by the — 1 functions COj{r). We comment 
that our ansatz ( fTOl i is by no means the only possible choice in su{N) EYM. Tech- 
niques for finding all spherically symmetric su{N) gauge potentials can be found in 
where all irreducible models are explicitly listed for < 6. 
With the ansatz ( fTOl i. there are N — I non-trivial Yang-Mills equations for the 
— 1 functions coj: 

( 2Ar'\ 
+ [ 2m - 2r^p0 — j co'j + WjCOj = (12) 

for y = 1 , . . . — 1, where a prime ' denotes d/dr, 
1 

47 



2 

and too = con = 0. The Einstein equations take the form 



(13) 
(14) 



, S' 2G 

m' = lxG + r^Pe, T = " ' ^l^) 
S r 

where 

G=l^cof. (16) 

j=i 

Altogether, then, we have + 1 ordinary differential equations for the A^ + 1 un- 
known functions m{r), S{r) and COj{r). The field equations ( 112115b are invariant 
under the transformation 

(Ojir) ^ -COj{r) (17) 
for each j independently, and also under the substitution; 

j^N-j. (18) 

We are interested in black hole solutions of the field equations (112115b . We as- 
sume there is a regular, non-extremal, black hole event horizon at r = r/,, where (r) 
has a single zero. This fixes the value of m(r/,) to be: 

2m(r,)=r,,-^. (19) 

However, the field equations ( 112115b are singular at the black hole event horizon 
r = r/j and at infinity r ^ We therefore need to impose boundary conditions 
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on the field variables m(r), S{r) and (Oj{r) at these singular points. When the cos- 
mological constant A is zero, local existence of solutions of the field equations in 
neighbourhoods of these singular points has been rigorously proved 199111241 . This 
proof can be extended to the case when the cosmological constant is negative ISl fTTI . 

We assume that the field variables (Oj{r), m{r) and S{r) have regular Taylor series 
expansions about r = r/,: 

m{r) = m{rh) + m'{rh) {r-rh)+0{r~ rhf ; 
(Oj{r) = (Oj{rh) + a]{ri,){r-rh) + 0{r-rhf; 
S{r)=S{rh)+S'{rh){r-rh) + 0{r-rh). (20) 

Setting /i(r/,) = in the Yang -Mills equations (fT2] | fixes the derivatives of the gauge 
field functions at the horizon: 

coKoO- "^-^"'^"-^"-^ (21) 

2m{rh)-2rlpe{rh) 5-^ 

Therefore the expansions ( l20b are determined by the 1 quantities a)j(r/,), r/,, 
5(r/,) for fixed cosmological constant A. For the event horizon to be non-extremal, 
it must be the case that 

2m' (r,,) = Irlpe (r,,) < 1 - Arf , (22) 

which weakly constrains the possible values of the gauge field functions (Oj{ri,) at 
the event horizon. Since the field equations (112115b are invariant under the transfor- 
mation (fTTl i. we may consider (Oj{rh) > without loss of generality. 

At infinity, we require that the field variables COj{r), m{r) and S{r) converge to 
constant values as r ^ 00, and have regular Taylor series expansions in ' near 
infinity: 

m{r)=M + 0{r-^)- 5(r) = 1 + O (r^^) ; C0j{r) = C0j,^ + 0{r-^) . (23) 

If the space-time is asymptotically flat, with A = 0, then the values of (Oj^co are 
constrained to be 

(Oj,^ = ±^j{N-j). (24) 

This condition means that the asymptotically flat black holes have no magnetic 
charge at infinity, or, in other words, these solutions have no global magnetic charge. 
Therefore, at infinity, they are indistinguishable from Schwarzschild black holes. 
However, if the cosmological constant is non-zero, so that the geometry approaches 
(a)dS at infinity, then there are no a priori constraints on the values of coj^oo. In 
general, therefore, the (a)dS black holes will be magnetically charged. It should be 
noted that the boundary conditions in the case when the cosmological constant A 
is positive are more complex, as there is a cosmological horizon between the event 
horizon and infinity. 
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Although the field equations (112115b are highly non-linear and rather complicated, 
they do have some trivial solutions which can easily be written down: 

Schwarzschild(-(a)dS) Setting 



COjir) EE ±V7(A^ (25) 

for all j gives the Schwarzschild(-(a)dS) black hole with 

m{r) = M = constant (26) 

We note that, by setting M = 0, pure Minkowski (A = 0) or (a)dS (A ^ 0) space 
is also a solution. 
Reissner-Nordstrom( -( a )dS) Setting 

COj{r)=0 (27) 

for all j gives the Reissner-Nordstrom(-(a)dS) black hole with metric function 

2M Q2 ^2 

M'-Hl + ^-^' (28) 

r 3 

where the magnetic charge Q is fixed by 

Q"^lNiN+l)iN-l). (29) 
o 

Embedded su{2) solutions For our later numerical and analytic work, an addi- 
tional special class of solutions turns out to be extremely useful. We begin by 
setting 

C0j{r)^±^j{N-j)(0{r) yj=l,...,N-l, (30) 

then follow ||99l and define 



XM=^^N{N-l)iN+l), (31) 

and then rescale the field variables as follows: 

R = X^\; A=A^A; m(i?) = 'm(r); 
S{R) = Sir); a){R) = Co{r). (32) 



Note that we rescale the cosmological constant A (this is not necessary in 
as there A = 0). The field equations satisfied by m{R), S{R) and &{R) are then 
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dm 
'dR 

1 dS _ IG 
SdR ~ ^"r"' 



HG + R^pg; 



O^R^^i 



d~cb 
dR^ 



2m — 2R^pq — 



2AR^ 



(33) 



where we now have 



and 



, 2m AR^ 

^ = '--R-—^ 



daV 

Ir) 



(34) 



(35) 



The equations {33[ are precisely the su(2) EYM field equations. Furthermore, 
the boundary conditions ( I20I23I ) also reduce to those for the su(2) case. 



2.3 Dyonic field equations 



As will be discussed in section 14.3.11 if either > 2 or we have a negative cos- 
mological constant A, then we do not need to restrict ourselves to considering only 
purely magnetic equilibrium gauge potentials. If the electric part of the gauge po- 
tential (|6]l, £/, is non-zero, there is still sufficient gauge freedom to set = in (|6]l 
||98l . Then, provided none of the cOj vanish identically, one of the Yang-Mills equa- 
tions again tells us that all the yj are identically zero. Following t98i it is convenient 
to define new real variables cx,j{r) by 



(36) 



so that the matrix £/ is automatically purely imaginary, diagonal and traceless. In 
this case the Yang-Mills equations ( fT2] i now take the form ||98ll 



„ / q 2Ar^\ , u 9 

lico'l + 2m - 2r^pe -— Oj'j + Wjdij + ^ajcoj 



0, 



and there are additional Yang-Mills equations for the a/, namely 



(37) 



(38) 



The Einstein equations retain the form ( fTsT i but the quantities pg ( flj] ) and G ( fTSI l 
now read ll98l 
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Pe 



4r4 ^ 



CO 



-N-l+2j) + 



.7=1 
N-l 

G= 

7=1 



-ajcoj]. 



(39) 



2.4 Perturbation equations 

We are also interested in the stability of the static, equilibrium solutions. For sim- 
plicity, we consider only linear, spherically symmetric perturbations of the purely 
magnetic solutions. We return to the general gauge potential of the form (|6]l, and 
the metric Q, where now all functions depend on time t as well as r. There is still 
sufficient gauge freedom to enable us to set = 0. This choice of gauge is partic- 
ularly useful as then we shall shortly see that the perturbation equations decouple 
into two sectors, the "gravitational" and "sphaleronic" sectors fllOll . We consider 
perturbations about the equilibrium solutions of the form 

(Oj (f , r) = (Oj {r) + d COj (r , r) , (40) 

where (Oj{r) are the equilibrium functions and 5cOj{t,r) are the Unear perturbations. 
There are similar perturbations for the other equilibrium quantities m and S, and in 
addition we have the perturbations dyj{t,r) and 5fij{t,r), the latter being the entries 
along the diagonal of the matrix 3§ ©I 

^ = Diag(/5j3i,...,/5j3w). (41) 

Note that the 5)3/ are not independent because the matrix ^ is traceless, so 

5/31 + .. . + 5/3^ = 0, (42) 

but it simplifies the derivation of the perturbation equations to retain all the 5)3/ for 
the moment. We ignore all terms involving squares or higher powers of the pertur- 
bations. The full derivation of the perturbation equations is highly involved and the 
details will be presented elsewhere IfTTI . Instead here we summarize the key features 
of the perturbation equations. As usual, we will employ the "tortoise" co-ordinate 
r*, defined by: 

^ = ^' (43) 
dr jlS 

where and S are the equilibrium metric functions. 
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The sphaleronic sector consists of the 2N — 1 perturbations 8^j, j = 1,. . . ,N and 
Syj, j — i, - ■ ■ ,N — 1. We define new variables 80] by 



5(Pj = cOjSyj. 



(44) 



The perturbation equations for the sphaleronic sector arise solely from the Yang- 
Mills equations, and comprise: 



[0)2 (5/3,+, - 5/3,) - co]_, (5/3, - 5i3,_i)] ; 



(45) 



dl {8<Pj) ^ {dl CO,) 50, + ^Scojd,, (5^, ~ 5^,+i) 
+ [ju idr,S) coj + idr,^)Sc0j + 2^iS{dr,C0j)] (5/3, - 5j3,+i) ; (46) 
together with the Gauss constraint 



= dr, 5j3, 



5$j + - [(Oj50j + a),_i5*,_i] , (47) 



where a dot denotes d/dt. It is important to note that the cosmological constant A 
only appears in these equations through the metric function jj. (|5]), and therefore the 
perturbation equations ( 145 146 1 ) and the Gauss constraint ( l47b have exactly the same 
form as derived in (47| for arbitrary gauge groups in asymptotically flat space. 



2.4.2 Gravitational sector 

The gravitational sector consists of the perturbations of the metric functions SjX and 
55 as well as the perturbations of the remaining gauge field functions 5a),-. Both 
the Einstein equations and the remaining Yang-Mills equations are involved in this 
sector. For an arbitrary gauge group and asymptotically flat space, the perturbation 
equations in this sector have been considered in BtI . In asymptotically adS, we 
also find that the metric perturbations can be eliminated to give a set of equations 
governing the perturbations 5o,, which can be written in matrix form 

M = dl (M) + .-^gM, (48) 
where 5w = {5coi 5cOn- i )^ and the (A^ — 1 ) x (A^ — 1 ) matrix has entries 
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uS^ 4 85" 

^GJ.j = V [Wj - 20)2] + —Y {dr,COjf + -WjCOj idr,COj) ; 

r iXof r 

uS^ 4 

-^GjJ+l = —0}jO}j+i + —r{dr,,0}j){dr,(Oj+l) 

ss 

+ [Wj(Oj {dr,(Oj+i) + Wj+i(Oj+i {dr,(Oi)] ; 

4 85 
-^GjM - — r(a,-,a),)(<9,-,a)t.) + — [W,«,(a,^a)A.) + WA.coA.(a,,a),)]; (49) 

where /: 7^ J , J + 1, and Y is given in terms of the equilibrium metric functions ji and 
S as follows; 

r=-dr,IX + ];dr,S+^. (50) 



3 Asymptotically flat/de Sitter solutions for su(A^) EYM 

We now turn to black hole solutions of the EYM field equations, beginning by briefly 
reviewing some of the key features of solutions in asymptotically flat or asymptoti- 
cally de Sitter space. 



3.1 Asymptotically flat, spherically symmetric su(2) solutions 

Apart from the trivial solutions given above ( 125127b . the first black hole solutions of 
the EYM field equations were found by Yasskin 01811 . and correspond to embed- 
ding the Reissner-Nordstrom electromagnetic gauge field into a higher-dimensional 
gauge group. The metric of these solutions is still Reissner-Nordstrom. Yasskin con- 
jectured that his solutions were the only ones possible. This conjecture was only 
shown to be false twenty-five years later lfT9l 1 1001 1 1671 1 1681 . That the discovery of 
hairy black holes in su(2) EYM took so long may be attributed to the conjecture that 
there were no soliton solutions in this model. This conjecture is based on the fact 
that there are no solitons in pure gravity (see, for example, 1178111031 ); no solitons in 
Einstein-Maxwell theory 177 1 , no pure YM solitons in flat space- time ||53ll56l and 
no EYM solitons in three space-time dimensions ||57l . However, once Bartnik and 
McKinnon [[7] had discovered non-trivial EYM solitons in four-dimensional space- 
time, Yasskin's no-hair conjecture for EYM theory was quickly shown to be false 

CD. 

For 5u(2) EYM, it has been shown ll23l l62l l67l that non-trivial solutions (i.e. 
solutions in which the gauge field is not essentially Abelian) must have a purely 
magnetic gauge potential, which is described by a single gauge field function (o{r) 
(fTOb . Note that the ansatz ( fTOl i for su(2) is not the same as the Witten ansatz II178I 
which was used in the original papers QIH], but it gives equivalent field equations. 
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In this case the su(2) EYM equations have the form 



1 dS 
'sTr 



dm 
"57 









(1-0)^)^] dO) 



+ [l-a)2]a). (51) 



r dr 



It is the highly non-linear nature of these equations which allows for non-trivial 
soliton and hairy black hole solutions, which may be thought of heuristically as 
arising from a balancing of the gravitational and gauge field interactions (see 1821 
for a recent discussion). The non-Unear nature of the equations also means, however, 
that (apart from the solutions for the Yang-Mills field on a fixed Schwarzschild 
metric ||28l[34| ) solutions can only be found numerically. 

The numerical work in 171 |T9l 1 1001 1 1671 [1681 found discrete famiUes of solutions 
II155I , indexed by the event horizon radius r/, (with r/, = for solitons) and «, the 
number of zeros of the single gauge field function ©, each pair identifying a 

solution of the field equations. A key feature of the solutions is that n > 0, so that 
the gauge field function must have at least one zero (or "node"). Later analytic work 
129 , 148 149, 150J rigorously proved these numerical features. The black holes are 
"hairy" in the sense that they have no magnetic charge I23ll62ll67l . and are therefore 
indistinguishable at infinity from a standard Schwarzschild black hole. However, the 
"hair", that is, the non-trivial structure in the matter fields, extends some way out 
from the event horizon, leading to the "no-short-hair" conjecture 11211 . 

Although initially controversial l20l l24l 11511 11721 . rapidly it was accepted that 
both the soliton 01531 and black hole solutions II 1541 are unstable. This instability is 
not unexpected if we consider the solutions as arising from a balancing of the gauge 
field and gravitational interactions. Studies of the non-linear stability of the solutions 
[182, 183] reveal that the gauge field "hair" either radiates away to infinity or falls 
down the black hole event horizon, leaving, as the end-point, a bald Schwarzschild 
black hole. Due to this instability, the black holes, while they violate the "letter" of 
the no-hair conjecture, may be thought of as not contradicting its "spirit", and one 
might be led to conjecture that all stable black holes are fixed by their mass, angular 
momentum, and conserved charges. 

Originally these hairy black holes were shown to be unstable using numerical 
techniques 01541 but the instability can also be shown analytically l68l 11691 . In 
the su(2) case, the perturbation equations ( I45I46I48I I simplify considerably. The 
sphaleronic sector reduces to a single equation (see section 14.21 below for further 
details) 



r (0' 




(52) 
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while, on eliminating the metric perturbations, the gravitational sector also has just 
one equation: 



The instability has been compared to that of the flat-space Yang-Mills sphaleron 
II169I , which has a single unstable mode. The situation is slightly more complicated 
here, due to the two sectors of perturbations. The sphaleronic sector certainly, as its 
name suggests, mimics the perturbations of the flat-space sphaleron. It can be shown 
II 16611 that the number of instabilities in the sphaleronic sector equals «, the number 
of zeros of the gauge field function co. The same is true in the gravitational sector, 
as conjectured in 01011 and can be shown using catastrophe theory, by considering 
the more general EYM-Higgs solutions HI 141 . The above concerns only spherically 
symmetric perturbations. It is known that the flat-space sphaleron has instabilities 
only in the spherically symmetric sector [4J. Extending this to the su(2) EYM black 
holes requires complicated analysis II142I . using a curvature-based formalism de- 
veloped in l43l[T43l[T44l . 

Using the isolated horizons formalism, these "hairy" black holes can be inter- 
preted as bound states of ordinary black holes with the Bartnik-MacKinnon solitons 
m |54] |55] . In particular, the sohton masses are given in terms of the masses of the 
corresponding black holes ll55l . and the instability of the coloured black holes arises 
naturally from the instability of the corresponding solitons 131 1541 . 

Since these initial discoveries a plethora of new, asymptotically flat, hairy black 
hole solutions have been found in Einstein- Yang-Mills theory and its variants (see 
II 1 701 for a review of those solutions discovered prior to 1999). Most of these are, 
indeed, unstable. However, there are notable exceptions, including: (a) the Skyrme 
black hole ll22l l58l [80l [HTl where the existence of an integer-valued topological 
winding number renders the solutions stable, (b) Einstein- Yang-Mills-Higgs black 
holes in the limit of infinitely strong coupling of the Higgs field H] and (c) a par- 
ticular branch of Einstein-non- Abelian-Proca black holes EH [l09l [l58l [BS] [iMl • 
We will not consider additional matter fields further in this article. 



3.2 Non-spherically symmetric, asymptotically flat su(2) solutions 

One of the surprising aspects of the failure of black hole uniqueness in EYM is 
that almost every step in the uniqueness theorem in Einstein-Maxwell theory has 
a counter-example in EYM (see 1791 for detailed discussions on this topic, and 
ll45l 11271 11521 11551 11561 for examples of some results from Einstein-Maxwell the- 
ory which do generalize). An important example of this is Israel's theorem 18611871 . 
which states that the geometry outside the event horizon of a static black hole must 



-5(0 ^ -dj^ (5(b) 



(53) 
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be spherically symmetric. This is not true in EYM: there are static black hole so- 
lutions which are not spherically symmetric but only axisymmetric ||90l (in more 
general matter models, static black holes do not necessarily possess any symmetries 
at all 11371 11381 ). These solutions are found numerically by writing the metric in 
isotropic co-ordinates: 

ds^ = -/ r, e)dt^ + -jT^dr^ + ;/ ' dO^ + ^' d^ , (54) 

/(r,0) /(r,0) /(r,0) 

and using the following ansatz for the su(2) gauge field II136I : 

A = 1 {t; [//i (r, 0) t/r + (1 - //.(r, 0)) rdO] 

-p [T!^Hi{r, e) + T'g{l- //4(r, 0))] rsin Odij)}, (55) 

where 

Tj.' = T^(sin0cos/90,sin0sinp0,cos0) , 

Tg = (cos cos p(p, COS sinp^, — sin 0) , 

T^' = T,(-sin/70,cosp^,O), (56) 

with 

T= (-r„Ty,T2), (57) 

where T^, Ty, T- are the usual generators of su(2). Here, p is a winding number, with 
p = I corresponding to spherically symmetric solutions (with the gauge potential 
written in a different form to that we have used in (fTOli). Substituting the ansatz into 
the field equations gives a complicated set of partial differential equations, solutions 
of which are exhibited in 1901 . Static, axisymmetric soliton solutions also exist ll65l 
[85] ED. 

It is less surprising that rotating black holes also exist in this model 1921 193]| . 
generalizing the Kerr-Newman metric (as predicted in 11551 ). These solutions are 
indexed by the winding number p ( |56] | and a node number n. They carry no magnetic 
charge, but all have non-zero electric charge II1551 11561 . The question of whether 
there are rotating solitons in pure su(2) EYM has yet to be conclusively settled, 
however. Rotating soliton solutions have been found in EYM-Higgs theory II126I . 
but not in pure EYM theory. Although rotating solitons are predicted perturbatively 
l44l . the consensus in the literature is now that it seems unhkely that rotating soliton 
solutions do exist ifTTl . 



3.3 Asymptotically flat su{N) solutions 



We shall next consider generalizations of the su(2) YM gauge group. The simplest 
such generalization is to consider 5u(A^) EYM. The results of ll62l l67l do not ex- 
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tend to this larger gauge group, and it is possible to have solutions with electric 
charge ||69l . which correspond to a superposition of electrically charged Reissner- 
Nordstrom and the su(2) EYM black holes. Numerical solutions of the field equa- 
tions have been found in the following papers: 1 69 94 95 96J. As increases, the 
possible structures of the gauge field potential (|6]l become ever more complicated. 
A method for computing all spherically symmetric su{N) gauge field potentials is 
given in |6j, where all the irreducible possibilities are enumerated for N < 6. As 
in the su(2) case, black hole solutions are found at discrete points in the parameter 
space {C0jin,),j = I... N-1}. 

There is comparatively little analytic work for more general gauge groups. Local 
existence of solutions of the field equations ( I12I15I ) near the black hole event hori- 
zon and at infinity has been proven for gauge group su(A^) 1991 . and subsequently 
extended to arbitrary compact gauge group II123II124I . The existence of non-trivial 
black hole solutions to the field equations has been proven rigorously only in the 
su(3) case II139II1401 . although there are arguments that hairy black hole solutions 
exist for all II115I . In the su(3) case, Ruan II 1391 11401 has proved that there are 
infinitely many hairy black hole solutions, indexed by the numbers of zeros («i , «2) 
respectively of the two gauge field functions (0)1,0)2). Furthermore, provided that 
the radius of the event horizon is sufficiently large, there is a black hole solution 
for any combination of («i,«2)- The global properties of the solutions for arbitrary 
compact gauge group are studied in |125|. However, it will come as no surprise 
to learn that all these solutions, in asymptotically flat space, and for any compact 
gauge group, are unstable 14611471 . To show instability it is sufficient to find a single 
unstable mode, and therefore the work in l46l l47l studies the simpler, sphaleronic 
sector of perturbations (see section l2".4.1l i. 



3.4 Asymptotically de Sitter su(2) EYM solutions 

Another natural generalization of asymptotically flat su(2) EYM is the inclusion of 
a non-zero cosmological constant A . When the cosmological constant is positive, 
soliton 11711 and black hole 11631 su(2) EYM solutions have been found (other 
numerical solutions are presented in BTl II 181 ). These solutions possess a cosmo- 
logical horizon and approach de Sitter space at infinity (for a complete classification 
of the possible space-time structures, see lISOl ). The phase space of solutions is again 
discrete, and the single gauge field function O) must have at least one zero. Unsur- 
prisingly, these solutions again turn out to be unstable l42l |63l [1631 . Given this 
instability, the asymptotically de Sitter solutions have received rather less attention 
in the literature, but some analytic work can be found in ||104||105lll06| . 
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4 Asymptotically anti-de Sitter solutions for su(2) EYM 



We now turn to the main focus of this article: asymptotically anti-de Sitter solutions. 
We begin by reviewing some of the properties of black holes in su(2) EYM. 



4.1 Spherically symmetric, asymptotically adS, su(2) EYM 
solutions 



studied in II174II . and subsequently in 



Black hole solutions of su(2) EYM with a negative cosmological constant were first 

The field equations now take the form 

2^2. 



dm 
'd7 

Ids 

S dr 
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CO 



2 fdcoy 

7 \~d7) ' 



= 
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2Ar^ 



1 - CO 



2^21 



dco 
dr 



(58) 



The inclusion of a negative cosmological constant means that boundary conditions 
at infinity (l23T l are considerably less stringent than in the asymptotically flat case; it 
is therefore unsurprising that it is easier to find solutions in asymptotically adS. 

The space of solutions in adS is very different to that in asymptotically flat space. 
Instead of finding solutions at discrete values of fi)(r/,), instead solutions exist in 
continuous, open intervals. Furthermore, for sufficiently large |A |, we now find so- 
lutions in which the single gauge field function co{r) has no zeros. A typical exam- 
ple of such a solution is shown in figure [T] further examples can be found in Iil74] . 
These properties of the space of solutions of the equations (|58] | are proved in II174L 

We now examine the structure of the space of solutions, more details of which 
can be found in ||8ll9l ll74l . There are three parameters describing the solutions, r/,, 
A and Co{ri,). In order to plot two-dimensional figures, we fix either r/, or A and 
vary the other two quantities. For su(2) black holes, the constraint (|22] ) on the value 
of the gauge field function at the event horizon reads 



{co{r,f-iy<rl{l-Arl) 



(59) 



Whether we are varying r/, or A, we perform a scan over all values of O)/, which 
satisfy ( |59] ). Firstly, we show in figure|2]the space of black hole solutions for fixed 
A = —0.01 and varying event horizon radius r;,. The outermost curves in figure|2] 
are where the inequality ( |59] l is saturated. Immediately inside these curves we have 
a shaded region, which represents values of (r/,, a)(r/,)) for which the constraint 
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0.0 -I 1 1 < < 1 

10° 10^ 10^ 10' 10" r 10= 

Fig. 1 An example of an su(2) EYM black hole in adS in which the gauge field function ff)(r) has 
no zeros. Here, A = — 1, r;, = 1 and (B(r;,) = 0.7. 




Fig. 2 The space of 5u(2) black hole solutions when A = —0.01, for varying ri,. The shaded region 
indicates values of the gauge field function ft)(r/,) at the event horizon for which the constraint \59) 
is satisfied, but for which we find no well-behaved black hole solution. The number of zeros n of 
the gauge field function ft) are indicated in those regions of the phase space where we find black 
hole solutions. Elsewhere on the diagram, the constraint (59) is not satisfied. Between the region 
where n = 2 and the shaded region we find black hole solutions with n = 3, 4 and 5, but these 
regions are too small to indicate on the graph. Taken from (9|. 
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( |59] ) is satisfied, but for which we are unable to find black hole solutions which 
remain regular all the way out to infinity. Where we do find solutions, we indicate in 
figure|2]the number of zeros of the gauge field function ft) (r) . The solution for which 
ft)(r/,) = 1 is simply the Schwarzschild-adS black hole, while that for ftj(r/,) = is 
the magnetically charged Reissner-Nordstrom-adS black hole (see section |272] |. As 
r/, — > 0, the constraint (l59T l implies that ftj(r/,) 1, as can be seen in figure |2] The 
black hole solutions become soli tons in this limit. However, for this value of A, 
there are different soli ton solutions, with ft) having different numbers of zeros lISTI . 
a feature which is not readily apparent from figure |2] We find similar behaviour on 
varying r/, for different values of A . 

If we now fix the event horizon radius to be r/, = 1 and vary A, the solution 
space is shown in figure [3] with a close-up for smaller values of \A \ in figure ID 




n = 4 n = 3 log,o(-A) 

Fig. 3 Phase space of su(2) black holes with r/, = 1 and varying A. The shaded region indicates 
values of the gauge field function ft)(r/,) at the event horizon for which the constraint )59t is sat- 
isfied, but for which we find no well-behaved black hole solution. The number of zeros n of the 
gauge field function m are indicated in those regions of the phase space where we find black hole 
solutions. Elsewhere on the diagram, the constraint (59} is not satisfied. As well as the regions 
where n = 0, ... ,4 as marked on the diagram, we find a small region in the bottom left of the plot 
where n = 5. This region is too small to indicate on the current figure, but can be seen in figure|4] 
Taken from (9). 

Again, in figures |3] and |4] we have shaded those regions where the constraint ( |59] l 
is satisfied, but no regular black hole solutions could be found. Where we do find 
solutions, the number of zeros of the gauge field function co{r) is indicated in the 
figures. As A 0, the phase space breaks up into discrete points, which correspond 
to the asymptotically flat 'colored' su(2) black holes described in section [STI 1191 . 
For sufficiently large |A |, we find solutions in which the gauge field function has no 
zeros. 
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Fig. 4 Close-up of the phase space of su(2) black holes with r;, = 1 and smaller values of A. In 
the bottom left of the plot there is a small region of solutions for which n = 7, but the region is too 
small to be visible. Taken from j9J. 



The spectrum of black hole solutions (that is, the relationship between the mass 
M and magnetic charge Q of the black holes) was first studied in |!26|. We plot in 
figure|5]the black hole mass versus magnetic charge for black holes with r/, = 1 and 
varying values of A (cf. figure 8 in ll26l ). For large values of |A|, there are only 




Fig. 5 Black hole mass M and magnetic charge Q for su(2) EYM black holes with r;, = 1 and 
varying A (cf. figure 8 in 1261 ). 
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nodeless solutions and the spectram is simple, with the black holes being uniquely 
specified by A, and Qm- As \A \ decreases, the spectrum becomes more compli- 
cated. For example, looking at the A = —0.1 curve in figure|5]we see that a branch 
structure emerges. The lower M curve for A — —0.1 consists of « = (nodeless) 
solutions, and extends from negative Q up to Q = 1. When Q — l,a branch of « = 1 
solutions appears, which have larger mass. As Q decreases along this branch of solu- 
tions, the mass M increases, until a bifurcation point is reached and a second branch 
of n = 1 solutions appears, with even larger mass, and with the charge increasing 
as M increases. For smaller values of \A\, we find ever more complicated spectra, 
which appear to become "fractal" as |A| ^ ll26l II 131 . In view of the catastro- 
phe theory analysis of other hairy black hole solutions 01581 11591 11621 . one might 
anticipate that the stability of the solutions changes at the points in the spectrum 
where two branches of solutions meet, but this has yet to be fully investigated in 
the literature (see 11311 for an in-depth stability analysis of the soliton solutions). We 
therefore next consider the stability of these black holes. 



4.2 Stability of the spherically symmetric solutions 



As discussed in section ISTI for the asymptotically flat su(2) EYM black holes, it 
has been shown that the number of instabilities is twice the number of zeros of 
the gauge field function 0}{r). Therefore, one might anticipate that at least some 
solutions when co{r) has no zeros could be stable. For the su(2) EYM case, the 
perturbation equations ( I45I46I48] | simplify considerably. In the sphaleronic sector, 
there is a single 8<P ( l44l i and two further perturbations 5j3i, 5)32, although these are 
not independent ( |42] |. so we may consider just 8v = SjSo — SPi. The sphaleronic 
sector perturbations equations ( I45I46I ) then reduce to 



(0 dv; 



50 = dl {50) - - {dl(o) 50-^S(odr, (5v) 



CO 



[ju {dr,S) (0 + {d,.^ii)S(0 - 2pLS (<9r, (o)] 5v; 



and the Gauss constraint (l47b is now 

= d,-, (5v; 



2^5 5^ 
r S 



5v + —co5<P. 



(60) 



(61) 



(62) 



By introducing a new variable ^ by (note our notation above is different from that 
used in lfT74l ) 



-5v. 



(63) 



the sphaleronic sector then reduces to a single equation Ml 741 
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^^(1 + 0)2) + — — 



while the gravitational sector (|48T l also has just one equation: 
- 5(a = (5to) 



.-3 



(64) 
(65) 



The sphaleronic sector equation ( |64] | is exactly the same as that in the asymptotically 
flat su(2) EYM case (|52] |. but the gravitational sector equation ( l53T l unsurprisingly 
is modified by the presence of non-zero A . Both equations ( l64b and ( |65] | have the 
standard Schrodinger form 



(66) 



with potential For the sphaleronic sector, when the gauge field function o(r) 
has no zeros, it is immediately clear that the potential is positive, so there are no 
instabilities in this sector (this result does not hold in the asymptotically flat case 
because the zeros of (nir) in that case mean that is not regular). The gravita- 
tional sector potential is more complex to analyze, but, for sufficiently large \A \ and 
a)(r/,) > 1/ \/3, it can be shown that the potential is positive and there are no insta- 
bilities in this sector either. Therefore there are at least some hairy black holes which 
are stable under linear, spherically symmetric, perturbations. It can further be proved 
that at least some of these solutions remain stable when non-spherically symmetric 
perturbations are considered 11451 1177II but the analysis is highly involved and so 
we do not attempt to summarize it here. 

It should be remarked that it is unhkely that all nodeless black hole solutions are 
stable, although this has not been investigated in the literature. An in-depth study 
of the corresponding solitonic solutions BTI has revealed that some soliton solu- 
tions for which a)(r) has no zeros, although they do not have any instabilities in the 
sphaleronic sector, do possess unstable modes in the gravitational sector. A scaling 
behaviour analysis of the solitonic solutions fS?! has shown that the stable soli- 
ton solutions can be approximated well by the stable solitons which exist on pure 
adS space. On the other hand, the unstable solitons are interpreted as the unstable 
Bartnik-MacKinnon solitons Q dressed with solitons on pure adS. 



4.3 Other asymptotically anti-de Sitter su(2) EYM solutions 
4.3.1 Dyonic solutions 

In asymptotically adS, it is no longer the case that the only genuinely non-Abelian 
solutions must have vanishing electric part in the gauge potential Q, so the results 
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of ||62ll67ll do not extend to non-asymptotically flat solutions. As well as the mag- 
netically charged solutions described above, dyonic black holes were discussed in 
||25ll26l . which we shall not consider further here. The stability of the dyonic solu- 
tions remains an open question as the perturbation equations do not decouple into 
two sectors in this case, making analysis difficult. 



4.3.2 Topological black holes 

As in Einstein-Maxwell theory, topological black hole solutions exist for su(2) 
EYM in adS lfT6l . The metric in this case reads 



where 



and 



ds^ ^ -ilS^ dt"- + pi-' dr^ + dB- + r Y {B)dip\ (67) 

{sin0 for k=\, 
9 for k = 0, (68) 
sinh 9 for A: = — 1 , 

2m(r) A-p- 

^=k . (69) 

r 3 



The ansatz for the purely magnetic gauge field potential is now lfT6l 

d\nf 



A = TxC0(r)d9' 



'VyCO{r) + T:- 



d9 



fi9)d<j). (70) 



When A ~0, only spherically symmetric solutions with k = I are possible, but for 
y\ < 0, solutions with both k = Q and k — ~l have been found ||T6| . All the solutions 
are nodeless, which can be easily proved from the field equations (TSl. It is found in 
|fT6ll that all the A: = solutions are stable under spherically symmetric perturbations 
in both the sphaleronic and gravitational sectors. The same is true for the k = — I 
solutions for which o) > 1 as r ^ oo ifTSl . 



4.3.3 Non-spherically symmetric solutions 

As in the asymptotically flat case, there are both soliton II 12811 and black hole 111351 
solutions which are static but not spherically symmetric, so that the metric and gauge 
potential take the form ( 154155b . Rotating black holes have also been found 111 121 . and 
there are also rotating dyonic soliton solutions 111301 . 
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5 Asymptotically anti-de Sitter solutions for su{N) EYM 

In the previous section we found that stable hairy black holes exist in su(2) EYM 
with a sufficiently large and negative cosmological constant. A natural question is 
therefore whether there are stable hairy black hole solutions of su(A^) EYM in adS, 
and we examine this question in this section. 



5.1 Spherically symmetric numerical solutions 

For any fixed A^, the field equations ( fT2] i and ( fTSl l can be solved numerically using 
standard techniques. We will outline briefly some of the key features of the black 
hole solutions for su(3) EYM. Details of the corresponding soliton solutions and 
the solution space for su(4) EYM can be found in ||9l|. 

For su(3) EYM, there are two gauge field functions coi (r) and (a2{r), and there- 
fore four parameters describing black hole solutions: r/,. A, (Bi(r/,) and (£)2{rh). Us- 
ing the symmetry of the field equations ( fTTI i. we set (0\ (rf,), fi>2('"/j) > without loss 
of generality. The constraint (|22] | on the values of the gauge field functions at the 
horizon becomes, in this case: 

[a)i(r,)2_2]'+ a>2(r,)2]'+ [2- a>2(r,,)2]' < 2rl (l ~A>j,) . (71) 

Two typical black hole solutions are shown in figures|6]and|2] The metric functions 




Fig. 6 Typical su(3) black hole solution, with r;, = 1, A = — 1, Ci)| (?■;,) = 1.2 and a)2(?';,) = 1.3. In 
this example, both gauge field functions have no zeros. Taken from |9). 



behave in a very similar way to the su(2) solutions, smoothly interpolating between 
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Fig. 7 Example of an su{3) black hole solution, with = 1, A = -0.0001, £Oi{r;,) = 1.184 and 
(Oiiri,) = 1.216. In this case, both gauge field functions have three zeros. Taken from (9j. 

their values at the horizon and at infinity. We note that S{r) in particular converges 
very rapidly to 1 as r ^ oo. In figure|6] we show an example of a black hole solution 
in which both gauge field functions have no zeros. We note that both gauge field 
functions are monotonic, however, one is monotonically increasing and the other 
monotonically decreasing. In our second example (figure|7ll both gauge field func- 
tions have three zeros. Although, in both our examples the two gauge field functions 
have the same number of zeros, we also find solutions where the two gauge field 
functions have different numbers of zeros (see figures [8] and |9]). 

We now examine the space of black hole solutions. Since we have four parame- 
ters, in order to produce two-dimensional figures, we need to fix two parameters in 
each case. We find that varying the event horizon radius produces similar behaviour 
to the su(2) case, so for the remainder of this section we fix r/, = 1 and consider the 
phase space for different, fixed values of A, scanning all values of a)i(r/,), Oii{rh) 
such that the constraint ( ItTI ) is satisfied. From the discussion in section |2l we have 
embedded su(2) black hole solutions when, from (l30l i: 



which occurs when (0\ (r/,) = COiirh)- 

In figureslSlflOlwe plot the phase space of solutions for fixed event horizon radius 
rh= I and varying cosmological constant A = —0.1, —1 and —5 respectively. In 
each of figures ISlfTOl we plot the dashed line (Oi [r^) = COiifi,), along which lie the 
embedded su(2) black holes. It is seen in all these figures that the solution space 
is symmetric about this line, as would be expected from the symmetry ( fTSl ) of the 
field equations. The solution space is found to be symmetric about the line ©i (r/,) = 
©2 (Oi ) not only in terms of where we find solutions, but also in terms of the numbers 



COi(r) = V2(0(r) = 0)2 (r) 



(72) 
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Fig. 9 Solution space for su(3) black holes with r;, = 1 and A = — 1. The shaded region indicates 
where the constraint illi is satisfied but we do not find black hole solutions. Outside the shaded 
region the constraint jlU does not hold. Where there are solutions, we have indicated the numbers 
of zeros of the gauge field functions within the different regions. For this value of A there is a large 
region in which both gauge field functions have no zeros. Taken from (9). 
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Fig. 10 Solution space for 5u(3) black holes with r;, = 1 and A = —5. It can be seen that for 
the vast majority of the phase space for which the constraint 17 It is satisfied, we have black hole 
solutions in which both gauge field functions have no zeros. Taken from (9). 



of zeros of the gauge field functions. To state this precisely, suppose that at the 
point ft)i(r/,) = fli, 0)2 (Oj) = fl2 we find a black hole solution in which (Bi(r) has 
«i zeros and a)2(r) has 112 zeros. Then, at the point (Bi(r) = 02, (02{r) = fli, we 
find a black hole solution in which CO\{r) has «2 zeros and a)i(r) has ni zeros. 
This is clearly seen in figures [8] and |9] and follows from the symmetry (fTSl l of the 
field equations. As we increase |A|, we find (see figures [SifTOli that the solution 
space expands as a proportion of the space of values of (Oi (ri,), fi)2(r/,) satisfying the 
constraint ( TTTI ). It can also be seen from figures [SlfTOl that the number of nodes of 
the gauge field functions decreases as \A \ increases, and that the space of solutions 
becomes simpler. For A = —0.1, there is a very small region of the solution space 
where both gauge field functions have no zeros. This region expands as we increase 
|A |, until for A = —5, both gauge field functions have no zeros for all the solutions 
we find. 

The solution space becomes progressively more complicated as increases, due 
to the increased number of parameters required to describe the solutions. However, 
the key feature described above is found; namely that for sufficiently large |A|, all 
the solutions we find are such that all the gauge field functions Wj have no zeros. 
These solutions are of particular interest since one might hope that at least some of 
them might be stable. 

As with the su(2) black holes we may consider the spectra of black hole solu- 
tions by plotting the relationship between the mass M and magnetic charge Q of 
the solutions (see figure |5] for the su(2) case). As may be expected, for higher 
the spectra are even more complicated than for su(2). In figure [TT] we plot some 
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of the possible values of M and Q for su(3) EYM black holes with A = —0.1 and 
r/, = 1 . In figure [TT] we have colour-coded the various possible numbers of zeros 



0, n2 = 

1, n2 = 




1.5 



2.5 



3.5 4 
Q 



Fig. 11 Black hole mass M versus magnetic charge Q for su(3) EYM black holes with ri, = I and 
A = —0.1. There are many different combinations of numbers of zeros of the gauge field functions 
(see figure [8]l, which are indicated by different colours. Here we have performed a scan over a 
grid of possible values of the gauge field functions at the event horizon, ft)i (r;,), ah{ri,), leading to 
discrete points in the spectrum. This is to enable the complicated structure of the spectrum to be 
seen. 



of the gauge field functions (cf. figure |8]l. We have used a discrete grid of initial 
values of the gauge field functions at the event horizon (ft)i (r/,), a)2(r/,)) and plot- 
ted discrete points so that at least some of the structure can be seen. In this case, 
because we have a four-parameter (A,r/,,a)i(r/,),a)2(r/i)) space of solutions of the 
field equations, even when A and r/, are fixed, we obtain two-dimensional regions 
in the {M,Q) plane, rather than curves as in the su(2) case. It can be seen from 
figure [TT] that the spectrum is very complicated, with the regions corresponding to 
different numbers of zeros of the gauge field functions overlapping. It is certainly 
the case that the black holes cannot be uniquely characterized by the four parameters 
(A,r/„M,e). 



5.2 Analytic work 

For any fixed value of A^, it is possible to examine the space of solutions numeri- 
cally. However, we would like to know whether there are solutions for all N, and, in 
particular, whether for all there are some solutions for which all the gauge field 



Classical Yang-Mills Black Hole Hair in anti-de Sitter Space 



29 



functions have no zeros, which we expect to be the case for sufficiently large \A\. 
Answering this question for general requires analytic rather than numerical work. 

In II174II . the existence of black hole solutions for which the gauge function co{r) 
had no zeros was proven analytically in the su(2) case. Since su(2) solutions can be 
embedded as su{N) solutions via dSOl l. we have automatically an analytic proof of 
the existence of nodeless 5u{N) EYM black holes in adS. However, these embedded 
solutions are "trivial" in the sense that they are described by just three parameters: 
rh, A and Co{rh). The question is therefore whether the existence of "non-trivial" 
(that is, genuinely su(A^)) solutions in which all the gauge field functions COj{r) have 
no zeros can be proven analytically. The answer to this question is affirmative, and 
involves a generalization to su(A^) of the continuity-type argument used in II174L 
The details are lengthy and will be presented elsewhere ifTTI . Here we simply outline 
the key steps in the proof. 

The main idea of the proof is sketched in figure [12] We wish to find black hole 



Black hole solution is 
regular everywhere if. 



Regular 
here.... 


....regular 
here 










..regular 

liroughout here 








r 



r = i-fj r = CO 



Fig. 12 Sketch of the main steps in the proof of the existence of non-trivial Bu(A') EYM black 
holes in adS for which all the gauge field functions have no zeros. We wish to find black hole 
solutions which are regular on the event horizon, regular everywhere outside the event horizon, 
and regular at infinity. We thank J. E. Baxter for providing this sketch. 

solutions which are regular on the event horizon, regular everywhere outside the 
event horizon, and regular at infinity. The proof proceeds via the following steps: 

1. We firstly prove (generalizing the analysis of ||99l to include A) that the field 
equations (I12I151 I and initial conditions at the event horizon ( l20b possess, locally 
in a neighborhood of the horizon, solutions which are analytic in r, r/,, A and 
the parameters (Oj{rh). As might be expected, the analysis of ||99| requires only 
minor modifications to include a negative cosmological constant. 

2. This enables us to prove that, in a sufficiently small neighborhood of any em- 
bedded su(2) solution in which (o{r) has no nodes, there exists (at least in a 
neighborhood of the event horizon) an s\x{N) solution in which all the (i)j{r) 
have no nodes. 
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3. Using the analyticity properties of the solutions of the field equations, we then 
show that these su(A^) solutions can be extended out to large ^ r/,, provided 
the initial parameters (Oj{ri^) are sufficiently close to those of an embedded su(2) 
solution in which (o{r) has no zeros. Furthermore, by analyticity, none of the 
COj{r) will have any zeros between the event horizon r/, and r^. 

4. The key part of the proof lies in then showing that these su(A^) solutions can 
be further extended out to r ^ oo and that they satisfy the boundary conditions 
( |23] l at infinity. This part of the analysis uses the properties of the Yang-Mills 
field equations (fT2] i in the asymptotically adS regime. As in the su(2) case 01741 . 
these have very different properties from the asymptotically flat case, and this 
makes it much easier to prove the existence of solutions. Furthermore, it can be 
shown that the gauge field functions COj{r) will have no zeros for r > ri. 

In summary, this process gives genuinely su(A') black hole solutions in which all 
the gauge field functions have no zeros, and which are characterized by the + 1 
parameters r/„ A and (Oj{rh). 



5.3 Stability analysis of the spherically symmetric solutions 

The remaining outstanding question is whether these new black holes, with poten- 
tially unbounded amounts of gauge field hair, are stable. We consider hnear, spher- 
ically symmetric perturbations only for simplicity. The analysis of II1451 11771 in 
the 5u(2) case revealed that, for sufficiently large |A|, stability under spherically 
symmetric perturbations continued to hold also for non-spherically symmetric per- 
turbations, and one might hope that a similar result will hold in the more complex 
su(A^) case. However, we leave this for future work. Even for spherically symmetric 
perturbations, the analysis is highly involved in the su{N) case and the details will 
be presented elsewhere ||8] |TT|. Here we briefly outline just the key features. The 
perturbation equations themselves can be found in section 1241 

5.3.1 Sphaleronic sector 

The sphaleronic sector consists of the perturbation equations (145146b together with 
the Gauss constraint ( |47] |. The analysis of this sector essentially follows that of 1471 
in the asymptotically flat case. We begin by defining yet more new variables, Sej, 
for; = l,...,A^by 

5Ej=r^5l5j, (73) 

then, after much algebra, the sphaleronic sector perturbation equations can be cast 
in the form 

-^=^^s^, (74) 
where the {2N — 1 ) -dimensional vector is defined by 
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'F={deu---,S£N,d<t>u---,d<t>N-i)- (75) 

and is a self-adjoint, second order, differential operator (involving derivatives 
with respect to r but not t), depending on the equilibrium functions COj{r), m{r) and 
S{r). The operator can be written as the sum of three parts. The first is of the 
form x^X for a particular first order differential operator % (whose precise form 
can be found in HKTT]) and is therefore manifestly positive and is regular if the 
gauge field functions O/ have no zeros. The second part vanishes when applied to a 
physical perturbation due to the Gauss constraint ( l47b . The third part is a matrix Y 
which does not contain any differential operators. It can be shown that the matrix 
is regular and positive definite provided the unperturbed gauge functions (Oj{r) 
have no zeros and satisfy the — 1 inequalities 




for all j ^\,. . .N —I, and all r>ri,. The inequalities ( |76] l define a non-empty sub- 
set of the parameter space. For example, we show in figure[T3]where the inequalities 
(|76] ) are satisfied for the gauge field functions at the event horizon, for the particular 
case of A = —10 and r/, = 1. From figure [13] we can see that there are some node- 




less solutions which satisfy the inequaUties ( |76] | at the event horizon. For any A^, it 
can also be proved analytically that, for sufficiently large |A|, there are non-trivial 
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su(A^) solutions, in a neighbourhood of some embedded su(2) solutions, such that 
the inequalities ( |76] l are satisfied at the event horizon. 

However, the requirements of (|76l ) are considerably stronger, as the inequalities 
have to be satisfied for all r>ri^. Our analytic work shows that, in fact, for any 
and sufficiently large |A |, there do exist solutions to the field equations for which the 
inequalities ( 176] ) are indeed satisfied for all r. This involves proving that for at least 
some solutions for which the gauge field function values at the event horizon lie 
within the region where the inequalities ( |76] | are satisfied, the gauge field functions 
remain within this open region. In figure[T4]we show an example of such a solution 
forsu(3)EYM. 



12 



S 10 



(02(r)^ against ia,(r)^ 

for the black hole solution 



(a^isf = 1+1/2 cf),(r)^ 



Fig. 14 An example of an 6u(3) solution for which the inequalities \16\ are satisfied for all r>ri,. 
In this example, A = — 10, r/, = 1 and the values of the gauge field functions at the event horizon 
are £0i (r/,) = 2, (Oiin,) = 1.95. Taken from (TO). 



5.3.2 Gravitational sector 

As might be expected, the gravitational sector perturbation equations (l48l l are more 
difficult to analyze than the sphaleronic sector perturbation equations. For stable 
solutions, we require the matrix ( l49b to be negative definite. For sufficiently 
large \A |, it can be shown that is indeed negative definite for embedded su(2) 
solutions, provided that (0^{r) > 1 for all r > r/, (the existence of such su(2) so- 
lutions is proved, for sufficiently large |A|, in 11741 ). As described in section 15^ 
above, our analytic work ensures the existence of genuinely su(A^) solutions in a 
sufficiently small neighborhood of these embedded su(2) solutions. These su{N) 
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solutions are such that the inequalities ( |76l ) are satisfied for all r > r/, (and therefore 
the solutions are stable under sphaleronic perturbations). The negativity of can 
then be extended to these genuinely su(A^) solutions using an analyticity argument, 
based on the nodal theorem of [[2] (see also 01771 for a similar argument for the non- 
spherically symmetric perturbations of the su(2) EYM black holes). The technical 
details of this argument will be presented elsewhere [n]. 

The conclusion of the work in this section is that there are at least some genuinely 
5u(A^) EYM black holes in adS, for sufficiently large |A|, for which all the gauge 
field functions (Hj have no zeros, and which are stable under spherically symmetric 
perturbations in both the sphaleronic and gravitational sectors. 



6 Summary and outlook 

In this review we have studied classical, hairy black hole solutions of su(A^) EYM 
theory, particularly spherically symmetric space-times and black holes in adS. We 
very briefly discussed some of the key aspects of the solutions in asymptotically flat 
space, which have been extensively reviewed in 117011 . Hairy black hole solutions 
exist for all A^, with — 1 gauge field degrees of freedom 11 15L however, all these 
solutions are unstable 1471 . Therefore, while these hairy black holes violate the "let- 
ter" of the no-hair conjecture (that is, their geometry is not completely fixed by 
global charges measurable at infinity), its "spirit" is maintained. In particular, stable 
equilibrium black holes are comparatively simple objects, described completely by 
just a few parameters. 

The main conclusion of this article is that this is not true in adS. The existence 
of stable hairy black holes in su(2) EYM 11741 did not really contradict the "spirit" 
of the no-hair conjecture, as only a single additional parameter was required to fix 
the geometry outside the event horizon. However, the recent work ifTOI which shows 
that there are stable hairy black holes in su(A^) EYM in adS for arbitrarily large 

changes the picture completely. For sufficiently large |A|, an infinite number of 
parameters are required in order to describe stable black holes. We might flippantly 
describe these as "furry" black holes, since they possess copious amounts of hair. 

What are the consequences for black hole physics in adS of these "furry" black 
holes? These need to be explored. Given the huge amount of interest in the adS/CFT 
correspondence in string theory 11 10II179]I180I . a natural question is how black hole 
hair in the bulk asymptotically adS space-time relates to the dual CFT. In particular, 
it has been suggested 1761 that there should be observables in the dual (deformed) 
CFT which are sensitive to the presence of black hole hair. Another example of this 
approach can be found in ifTOl , where an adS/CFT interpretation is given of some 
stable seven-dimensional black holes withso(5) gauge fields. We would expect that, 
in analogy with the su(2) case Il49ll50ll74l[84l [Till [1291 [T3T1 . there are solutions in 
some super-gravity theories with a gauge group containing an su{N) factor, which 
will need to be studied in the context of adS/CFT. There is evidence II 1611 that there 
are non-trivial black hole solutions of su(oo) EYM in adS, giving black holes not 
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just with unbounded amounts of hair, but infinite amounts of hair, at least in the 
hmit \A \ ao. It remains to be seen whether exact solutions of the 5u(oo) field 
equations can be found for finite A < 0, and whether any of these black holes are 
stable. If so, then their role in adS/CFT would be puzzling indeed. 

Due to space restrictions, there are many aspects of black holes in EYM which 
we have not been able to discuss. In particular, we have not mentioned the vast 
number of solutions which involve modifications of the EYM action ([TJ, including 
higher curvature terms (see, for example, ll88l[89l ) or the inclusion of dilaton (see, 
for example, 11331 ). Higgs (see, for example, lfT5lll07lll08ll ) or other modifications 
of the EYM action (see, for example. 111 191 11461 1147II ). Here we have also only 
studied four-dimensional space-times, while recent work has considered EYM in 
higher-dimensional space-times (see, for example, Il32l[33ll35ll36ll37ll38l[39ll40l 
ESI [ml [ml El and da for a review). 

The black hole solutions of EYM and its variants certainly exhibit an abundantly 
rich structure, and no doubt will have more surprises in store for us in the future. 
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